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SECTION  1 
INTRODUCTION 

1.1  Reason  for  Che  Study  and  Statement:  of  Problem  Conditions 

The  Ballistics  Research  Laboratory  (BRL)  has  developed  an  iterative 
method  to  determine  i.bital  parameters.  A  need  exists  to  have  an  initial 
approximation  correct  t.  approximately  10  per  cent.  This  study  has  been 
undertaken  by  the  Fhilco  Western  Development  Laboratories  to  determine  a 
method  for  finding  initial  values.  Errors  of  less  than  50  to  75  miles 
in  position,  and  1/2  to  1  mile/sec  in  velocity  are  specified. 

The  following  conditions  were  specified  for  the  system: 

a.  Method:  Polystation  Doppler 

b.  Ability:  Acquire  and  track  a  body  to  1000  miles  altitude. 

c.  Limitations: 

(1)  No  knowledge  of  vehicle  position  prior  to  acquisition 

(2)  No  assumption  about  the  nature  of  the  path 

(3)  No  cooperative  equipment  aboard  the  vehicle. 

1 . 2  Scope 

An  error  study  was  performed  to  determine  the  dependence  of  rms 
error  in  position  and  velocity  on  rms  error  in  the  observed  quantities. 

The  latter  were  assumed  equal  and  independent.  The  error  matrix  was 
calculated  for  a  number  of  circular  orbits  passing  through  points  in  a 
pre-assigned  grid. 

A  mathematical  model  was  constructed  for  data  reduction. 

Position-finding  consists  of  a  system  of  fourteen  linear  equations 
solved  for  approximate  values  of  the  position  coordinates.  These  are 
then  adjusted  by  the  standard  least  square  method  to  their  maximum- 
liklihood  values.  For  technique,  see  flow  chart  (Fig.  2a) 

For  velocity,  the  interval  immediately  following  the  vehicle's 
passage  through  the  grid  plane  is  considered.  For  this  interval,  readings 
to  two  stations  are  used,  in  conjunction  with  the  adjusted  position 
coordinates,  to  find  approximate  values  of  the  horizontal  and  vertical 
velocity  components.  These  are  then  adjusted  to  their  maximum- liklihood 
values  by  the  use  of  least  squares  technique  on  all  four  of  the  readings 
available  for  this  interval. 
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Test  cases  were  designed  for  the  mathematical  model  using  circular 
orbits  passing  through  points  in  a  grid  in  the  fan  beam's  central  plane. 
The  orbit  plane  inclination  was  varied  through  several  different  values 
(see  Section  48  and  Fig.  6) .  For  sequence  of  steps  see  flow  chart, 

?ig  3. 
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SECTION  2 

RESULTS  AND  CONCLUSIONS 

The  error  study  made  indicates  that  the  BRL  conditions  can  be  met 
Ky  the  proposed  system.  Errors  of  the  order  of  ten  wavelengths  can  be 
tolerated  at  300  me  The  error  study  which  was  performed  for  circular 
orbits,  resulted  in  data  which  is  presented  in  the  form  of  graphs,  some 
as  curves  and  some  as  straight-line  graphs 

To  permit  selection  of  optimum  beamwidth,  errors  were  studied  for 
values  trora  8°  to  30°  For  evaluation  of  the  selected  case,  charts  are 
presented  for  constant  beamwidths  In  order  to  aid  in  evaluating  a  system 
after  beamwidth  has  been  decided  upon,  graphs  of  error  propagation  factors 
were  plotted  holding  beamwidth  constant 
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SECTION  3 

RECOMMENDATIONS  FOR  FUTURE  WORK 

First  Ip  priority  for  further  work  would  be  continuation  of  the  test 

case  computation  lor  the  data  reduction  model,  (Sec  Section  6  ) 

Second,  as  a  leso  major  task,  a  more  realistic  error  calculation  for 

the  special  case  of  minimum-time- m-bcam  orbits  should  be  done 

Third,  it  would  be  desirable  to  vary  the  number  of  time  increments 

This  might  decrease  the  accuracy  of  the  approximate  values  determined, 

but  would  improve  the  least  squares  adjustments  There  is  reason  to 

believe  a  net  improvement  would  result 

Fourth,  a  review  of  the  conditions  imposed  on  the  system  reveals 

that  the  most  stringent  are  embodied  in  the  following  two  assumptions 

a  Passage  through  the  beam's  mid  plane  is  assumed  to  occur  at 
mid- time 

b  Time  of  entry  into  aid  exit  from  the  inner  beam  are  assumed  known 
In  the  present  system  it  is  the  coordinates  of  this  'piercing  point' 
for  which  the  system  is  solved  Unless  the  vehicle  is  so  cooperative  as 
to  maintain  a  velocity  configuration  and  reflection  characteristics  so  that 
it  actually  does  pierce  the  mid-plane  at  the  mid-time  between  first  and 
last  discernible  signals,  there  is  no  reason  to  suppose  that  this  assumption 
will  hold  in  an  operational  system  Of  course,  this  point  is  always  dis¬ 
cernible  if  a  split  beam  is  employed  However,  in  the  process  of  writing 
the  final  report,  it  was  determined  hat  both  of  these  conditions  could  be 
removed  for  a  fan  beam  as  well  and  thus  permit  the  solution  of  the  system 
for  any  point  on  the  path  of  the  vehicle  as  it  travels  through  the  beam 
Further,  it  is  possible  to  commence  and  cease  taking  readings  at  will, 
rather  than  having  to  assume  that  reading  commence  simultaneously  with 
entry  into  the  beam  and  cease  simultaneously  with  exit  from  the  beam  Un¬ 
certain  boundrv  knowledge  of  the  beam  renders  this  technique  questionable 
For  the  above  reasons,  it  is  strongly  suggested  that  the  scope  of  ike 
present  task  be  extended  to  include  a  thorough  investigation  of  the  above 
system. 
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SECTION  4 

SYSTEM  DESCRIPTION 


4  1  General 

The  Pcivstatlon  Doppler  Svsi em  has  the  ability  to  acquire  and  track 
a  sieving  vehicle  A  system  is  raid  to  acquire  if  it  can  determine  the 
position  of  a  vody  without  prior  knowledge  of  its  position  or  velocity. 

A  system  is  said  to  track  if  it  can  find  successive  points  on  a  vehicle's 
path  using  prior  position  information  A  Doppler  Polystation  System  as 
used  herein  refers  to  an  acquisition  and  tracking  system  of  more  than  one 
station  which  ttathers  only  information  regarding  diminution  or  augmentation 
of  received  cycles  due  to  motion  of  the  vehicle  relative  to  the  station 
complex 

This  particular  system  must  be  able  to  acquire  and  track  any  vehicle 
up  to  1000. miles  in  altitude.  The  region  of  coverage  is  300  miles  long 
and  1000  miles  high,  and  is  limited  laterally  by  the  sides  of  a  wedge- 
shaped  beam. 

4.2  Configuration 

The  system  is  composed  of  four  receivers  and  one  transmitter,  located 
along  an  arc  of  a  great  circle,  with  the  transmitter  at  one  end  of  the 
array  All  RF  signals  are  transmitted  from  the  ground,  reflected  from 
the  vehicle  and  received  at  four  ground  stations.  No  equipment  is  required 
ubouid  the  vehicle  A  computet  to  t educe  data  from  the  tracking  system 
is  the  only  ocher  equipment  required. 

4 . 3  Operation 

A  coninuous  signal  of  constant  frequency  is  transmitted  from  a 
ground  station.  This  signal  is  reflected  from  the  vehicle  and  received 
by  the  several  receivers.  The  time  during  which  signals  are  received  by 
all  four  receivers  is  the  period  of  observation  for  collecting  data  to 
be  reduced  for  position  and  velocity  information.  The  distance  from  the 
transmitter  to  the  vehicle  plus  the  return  to  the  receiver  is  defined  as 
the  range- sum. 

During  this  interval  a  certain  number  of  cycles  are  received.  This 
is  compared  with  the  number  of  cycles  transmitted  in  the  same 
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time  interval  Alternatively,  having  determined  che  total  time  of 

observation  -k  n.r,  d*--id-,  >t  into  a  number  of  sub-tntervals  and  determine 

the  cycle  count  difference  for  these  time  increments  Wo  shall  do  this  in 

the  following  way  Take  the  total  time  of  signal  reception  Divide  this 

into  fourteen  parts  Let  the  end  of  the  seventh  reading  be  the  reference 

point  r«l<  I’1*1  tnc  cycle  court  differences  received  in  intervals  of 

time  endir-g,  <  foi  rhe  fi  >t  half)  or  beginning  l  for  the  second  half)  at 

this  reference  time  The  so  will  be  denoted  hy  !}  .  They  will  represent 

cb  1 J 

the  range- sum  differences  between  the  1  time  and  the  reference  time 
measured  in  wavelengths;  tent  is,  »£ t j  is  the  increase  or  decrease  in  che 
range-sums  Ire  difference  betweer  the  transmitted  cycle  count  and 
received  cycle  count  of  the  i  time  -and  the  reference  time  (c  -  0)  is 
determined  by  using  a  common  time  reference  The  data  is  presented  to  a 
computer  which  solves  for  vehicle  position  and  velocity  The  time  at  the 
center  of  tin.  totar  tune  of  observation  is  detc. mined  and  position  and 
velocity  for  that  rime  are  determined  Knowing  this  point,  position 
coordinates  are  easily  determined  for  the  ocher  times 
4  4  System  Synthesis  Assumptions  and  Definitions 

The  operational  requirements  of  this  system  do  not  demand  any  assump¬ 
tions  about  the  -hspe  o£  the  earth  or  its  rotation  The  only  requirement 
is  chat  the  five  stations  be  co-planar  In  synthesizing  a  system  for 
discussion  in  rhis  report  a  spherical  earth  was  assumed  for  computational 
convenience  in  placing  the-  ground  stations  along  a  great  circle  In  the 
study  of  ttror  as  o  function  of  grid  position  lg-dop; ,  circular  orbits 
were  chosen  because  they  were  sufficiently  general  to  permit  reliable  con¬ 
clusions  to  be  drawn,  yet  not  so  complex  as  to  hamper  the  study  or  draw 
attention  from  the  true  purpose.  The  term  "cycle  count"  is  defined  as  the 
number  of  dopplc-r  cycles  counted  m  some  definite  time  selected  for  measure¬ 
ment,  as  opposed  to  "frequency"  which  usually  indicates  the  number  of 
cycles  in  one  second 

Position  is  referenced  to  a  right-handed  Cartesian  coord  ..rate  system, 
rotating  with  the  earth,  centered  at  the  transmitter,  and  oriented  as 
follows.  (1)  The  x-a^ls  passes  through  the  transmitter  and  the  last 
receiver,  positive  in  the  direction  of  the  receiver;  1 2 )  che  z-axis  lies 
in  the  plane  of  the  stations  and  passes  through  the  transmitter,  positive 
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"upward"  relative  to  the  earth;  (3)  the  y-axis  is  chosen  to  complete  a 
right-handed  Cartesian  system. 

Approximate  position  and  velocity  data  are  determined  by  the  solution 
of  certain  condition  equations.  It  is  assumed  that  these  values  sufficiently 
.oincide  with  the  true  values  to  warrant  a  first  order  approximation.  These 
-uadjusted  solutions  are  used  as  input  to  the  least  squares  adjustment  tech¬ 
nique  This  does  not  mean  that  this  is  an  iterative  technique  One,  and 
only  one,  adjustment  in  values  is  made. 

The  effects  of  both  special  and  general  relativity  were  studied.  In 
both  cases,  the  effects  were  found  to  be  non-existent.  Special  relativistic 
effects  are  cancelled  since  the  transmitter  and  receivers  do  not  move  relative 
to  each  other  The  possible  introduction  of  an  effect  due  to  accelerations 
of  the  vehicle  relative  to  the  station  complex  are  nullified  by  assuming  that 
the  energy  is  reflected  from  the  vehicle  in  zero  time. 

General  relativistic  effects  are  also  non-existent  In  the  proposed 
system  the  stations  are  all  located  on  the  surface  of  a  spherical  earth. 
Therefore,  all  stations  are  at  the  same  gravitational  potential.  While 
this  will  not  be  absolutely  true  for  any  operational  system,  nevertheless 
the  anticipated  variation  in  gravity  potential  at  different  station  locations 
would  be  beyond  the  precision  of  measurement  in  the  foreseeable  future, 
u  5  System  Synthesis 

Having  defined  the  coordinate  system,  the  coordinates  of  the  grouid 
stations  may  be  defined  as  follows.  (1)  Successive  ground  stations  are 
located  so  that  the  differences  of  their  x  coordinates  are  equal.  (2)  Since 
the  ground  stations  lie  along  a  great  circle  and  the  x-axis  is  defined  as  a 
cliuid  inteisecLiug  the  circle  at  the  transmitter  and  the  most  remote  re¬ 
ceiver,  the  z-coordinates  are  determined  by  knowing  the  radius  of  the  earth, 
see  Appendix  B  4  1  >'  )  The  y-coordinates  are,  of  course,  zero.  Thus, 

viewed  as  points  in  the  x- z  plane,  the  transmitter  is  located  at  (0,0)  ar.d 
the  four  receivers  at  R^(a,b),  R2(2a,c),  R^Oa.b),  R^(4a,0).  (See  Appendix 
B.4.1  for  derivation  of  expressions  for  b,c.)  (See  Figure  5.)  The  only 
condition  rec^uj red  b£  thi s  system  is  that  the  five  stations  be  co-planar 
For  computational  convenience  we  assume  they  are  located  along  a  great 
circle  on  the  surface  of  a  spherical  earth. 
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In  our  case,  che  x-coordinate  difference  referred  Co  in  (i)  was 
caken  at  200  miles.  The  resulting  values  of  z  were,  respectively,  for 
che  four  receiving  stations,  approximately  15,  20,  15,  0  miles.  Exact 
figures  are  derived  in  Appendix  fi  A  l.  The  more  accurate  figures  were 
used  in  all  computations  in  order  to  preserve  consistency  with  an  assumed 
spherical  earth  of  4000  mile  radius. 

The  direction  and  beamwidch  of  transmission  and  reception  create  a 
region  of  observation  which  is  wedge  shaped,  with  the  center  plane  of 
symmetry  being  che  sz  plane  The  edge  of  the  outermost  wedge  is  the  x-axis, 
the  inner  beam  edge  being  tangent  to  the  sphere  at  receiver  three  (See 
Fig  4  i  The  generating  angle  of  the  wedge  is  the  beamwidlh, 

4  6  Selection  ol  Optimum  Beamwidch 

To  pernii*-  the  selection  of  an  optimum  beamwidch,  sigma  matrices  for 
position  and  velocity  were  calculated  for  beamwidths  of  8°,  10°.  15°,  20°, 

OOO  0  ■ 

O  O  v  z  X  2 

25  ,  and  30  .  Charts  of  —  ,  —  ,  —  ,  and  —  are  presented  for 

°M  °M  °M  °M 

contours  of  constant  altitudes  for  various  combinations  of  horizontal 
displacement  xq  and  orbital-plane  angle.  (See  Fig  4  and  Figs  6  to  14.) 
o  is  the  standard  deviation  and  M  is  che  range-sum  difference. 

4  7  Error  Propagation  for  Fixed  Beamdwidth 

Having  selected  a  beamwidch,  Figs  15  to  20,  showing  the  behavior  of 
o  o 

X  2 

—  ,  —  ,  are  presented  as  functions  of  altitude,  holding  constant  x 
°M  0 

(horizontal  displacement)  and  0  Figures  21  to  28  are  also  presented 

a  a 
X  z 

showing  the  behavior  of  —  ,  —  ,  as  functions  of  horizontal  displacement 

°M  °M 

for  constant  altitude  and  0.  These  charts  are  for  constant  beamwidths 
and  are  presented  for  values  of  8°  and  15°.  By  super-imposing  on  a  light 
cable,  many  of  these  8°  to  15°  pairs  of  charts  may  be  compared  and  inter¬ 
polations  performed  to  find  a  values  for  intermediate  values  of  beamwidch 
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4 , 8  Statistical  Evaluation  of  the  System 

In  order  to  demonstrate  the  statistical  practicability  of  solving 
for  position  and  velocity  within  the  error  bounds  stipulated,  a  number 
of  representative  paths  were  synthetically  generated.  These  paths  pass 
through  the  beam  at  35  places,  intersecting  the  grid  plane  (the  xz  plane) 
at  angles  of  0,  15,  30,  45,  60,  75,  90  degrees  at  each  grid  point.*  In 
order  to  synthesize  a  situation  representative  of  operational  conditions, 
circular  orbits  were  selected  which  intersected  the  grid  plane  at  points 
defined  by  distances  of  50,  200,  350,  500,  and  650  statute  miles  along 
the  x-axis  and  distances  of  150,  300,  450,  600,  750,  900,  and  1050  statute 
miles  along  the  z-axis.  (See  Fig.  4  )  The  data  resulting  from  this 

study  are  presented  throughout  the  text.  The  equations  of  the  model 
are  included  in  Appendix  B. 


*  Those  for  90°,  of  course,  lie  in  the 
grid  plane. 
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SECTION  5 
DATA  REDUCTION 

Approximate  position  and  velocity  data  are  determined  by  the  solution 
of  certain  condition  equations-  It  is  assumed  that  these  values  coincide 
sufficiently  with  the  true  values  to  warrant  a  first  order  approximation. 
This  is  equivalent  to  saying  that  the  value  of  the  variance-covariance 
matrix  does  not  change  appreciably  in  the  region  containing  both  the 
adjusted  and  true  values. 

This  does  not  mean  that  the  mathematical  model  requires  an  initial 
point  or  prior  knowledge  of  the  vehicle's  "’ovement  It  means  that  the 
model  only  requires  that  the  (as  yet  unadjusted)  solutions  to  the  model's 
equations  be  sufficiently  free  from  arithmetic  error  propagation  to  warrant 
the  first  order  approximations , for  it  is  these  unadjusted  solutions  that 
are  used  as  input  to  the  leas t  squares  adjustment  technique. 

It  has  not  been  implied  that  this  is  an  .iterative  technique.  One, 
and  only  one,  adjustment  in  values  is  made.  If  the  unadjusted  value  is 
close  enough  to  the  correct  values  to  warrant  the  first  order  assumptions 
made,  the  adjusted  value  may  be  shown  to  be  the  maximum-likelihood  value 
in  a  least -squares  sense 

5  1  Wavelength 

In  this  system  wavelength,  a.  is  assumed  to  be  known,  since,  the 
transmitter  is  on  the  ground. 

5  2  Time  Duration  of  Cycle  Count 

The  time  of  entry  is  defined  as  the  first  point  in  the  beam  at  which 
the  vehicle  is  simultaneously  observed  by  all  four  receivers  We  shall 
assign  as  the  time  of  exit  from  the  beam,  a  point  which  is  symmetric  to 
the  entry  point  on  the  opposite  side  of  the  grid  plane  Zero  tin.-  will 
be  the  time  of  passage  through  the  grid  plane,  the  centermost  plane  of  the 
beam . 

5  •  3  Cycle  Count 

Information  is  received  in  the  form  of  a  plot  of  doppler  frequency 

in  cycles/sec  vs  time  in  seconds  at  each  of  the  four  receiver  locations 

for  the  full  period  of  time  (t.  to  t  )  that  the  vehicle  is  in  the  beam 

m  out 
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o£  antenna  2  This  may  then  be  integrated  tc  yield  total  cycles  ccunteu 
in  any  given  Interval 

Alternatively,  cycles  may  be  ccunted  directly,  thus  eliminating 
several  steps.  It  may  be  possible  to  utilise  the  computer  itself  as  a 
•.ounter  This  would  be  especially  attractive  since  many  modern  computers 
could  serve  as  a  self-contained  time  reference  while  accepting  and  pro¬ 
cessing  data  in  real  time 
5  A  Station  iccaticr. 

It  is  assumed  that  the  position  of  the  receiver  stations  1X^,2^)  relative 
to  cne  another  and  tc  the  transmitter  are  known  to  the  necessary  accuracy. 

As  described  earlier  (see  Section  A,  System  Description)  all  results  are 
relative  tc  a  c-crdinate  system  moving  with  the  stations 
5  5  Range- Sum  Differences  frem  Cycles  Counted 

The  observed  values  of  range-sum  differences  will  be  referred  to 

as  M. ,  or  M,  and  will  be  determined  as  follows: 
l)  k 


M 


c  (i-8)At  •  Qtj  \  i  -  1,  .  7 

c  (i-7)At  -  Q  X  i  =  8,  1A 

J  »  1,  •••  * 


where  c  Is  the  velocity  cf  propagation  in  miles  per  second  and  \  is 
the  wavelength  of  the  transmitted  frequency  in  miles 
When  written  with  a  single  subscript  =• 


for  k=j+4(i-l)  j  =  1,  L 
and  k  -  1 ,  ,56  i  =  1 ,  1A 


5  6  Pcsiticn  Determination 

561  Range- Sums  (Approximate)  From  Range-Sum  Differences 

Let  o  .  denote  the  distance  frem  the  transmitter  to  the  it*‘ 

1  •  th 

vehicle  pcsiticn  to  the  ^  station;  r.  ,  the  distance  from  the 

1C  th 

transmitter  tc  the  vehicle  position;  r  ^ ,  the  distance  from  the  i1"  vehicle 
position  tc  the  )rl1  receiver;  (X^,  Z^)  the  coordinates  of  the  jth 
ground  station  (note  that  the  transmitter  is  at  the  origin);  and 
(x^,  zr)  the  coordinates  cf  any  point  in  the  grid  plane 
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as  follows 


(3)  -2  r. 

10 

space 


°ij  =  rio  +  rij 
°ij  lio  “  rij 


Therefore 


r,  *  x  +  y.  *  z. 
io  i  i 


Let 


r.j2  =  (x^  x  ,2  <y.'  +  U.-Z.)2 

d.2«X*  **.2 
J  J  J 


Then 


.  2  ,  2 
(  o.  .  •  r.  j  =  (r.  .1 

lj  10  1J 


and  therefore 
2 


0. ,  •  2  r,  o  .  +  r.  2  =  r.  2  -2  X.  x.  -2  Z.z.  +  d  2 

lj  io  ij  io  io  J  l  J  l  j 


and  so 


2  X.x.  *  2  Z .i.  -  2  o. .  r.  =  a.  2  •  d.2 

J  i  J  r  iJ  ‘o  ‘J  J 


Four  equations  of  this  form  may  be  written  in  determinant  form 

to  ,  is  written  as  o.  for  convenience) 

°j  J 


The  determinant  is  zero  since  2z^  column  (4)  +  2  x^  *  column 
column  '2>  =  column  (1) 

Now  this  equation  is  valid  for  any  point  fx^,  y^,  z^,j  in 
When  this  is  written  for  the  i^1  vehicle  position, 


2 

(o1  +  Mu) 

0\  +  Mil 

xl 

Z1 

2 

(°2^i2>  • 

°2  +  Mi2 

x2 

Z2 

(03+Mi3»2  ‘ 

V 

p3  +  Mi3 

X3 

Z3 

<04  +  ‘ 

\2 

P4  +  Mia 

X4 

za 
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By  subtraction,  g.  =  f.  -f  =  0  which  may  be  written  in  the 
tern  ‘ 

h  *  l  %  “o  -  h  '  0 

P“1 

(These  will  be  used  as  condition  equations  for  the  least  squares  adjustment.) 


'l  = 

pl 

U5  =  pl 

u9  -  P, 

p2 

'2  = 

p2 

U6  =  p22 

„  2 

U10=  pl 

p3 

'3  = 

o3 

7  =  p3 

2 

UU=  P1 

p4 

4 

^4 

U8  =  p4 

U12=  P2 

P3 

The  solution  U^,  ...  of  these  fourteen  "pseudo"  linear 
equations  will  be  treated  as  approximations  Oj°,  ,  p3  ,  p^°  to  the  true  p.  . 

5.6.2  Range-Sums  to  Rectilinear  .Grid  Coordinates 

Using  three  (say  pk>  p.,  pm)  out  of  four  of  che  p^0,  ...  p^°, 
the  values  xq°,  z^0  and  incidentally  r^0,  may  be  found  by  solving 

2\  Xo°  +2Zk  z0°  V  rc°  -  (V2  •  <Pk°>2 

2:{1  Xo°  +2Z1  C  V  ro°  -  (dl>2  -  (Pi°)2 

2X  x  °  +2Z  z  °  -2p  °  r  °  =  (d  )2  -  (p  °)2 

mo  mo  Hmo  m  SKm 


For  example 

i  jk2  -(»k 


,2  ,  o  2  „  o 

dk  ~(pk>  Zk  r'k 


.  2  o2  o 

dl  '  (pl )  Z1  P1 

,  2  o2  ,  o 

d  -  (p  )  Z  p 

m  m  m  ^m 


X  Z 

m  m 


(It  is  not  suggested  Cramer's  rule  be  applied  in  actual  cases.  Modern 

matrix  inversion  techniques  almost  invariably  yield  more  reliable  results.) 
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The  "best"  set  C[cj  c^.  Pm  may  be  determined  by  evaluating 

4 

the  Jacobian,  J,  of  the  transformation  for  the  =  4  possible  combinations 

,0000 
Of  pL  c2  c3  PA 


5X0° 

az  0 

0 

0 

0 

c  1 

7*  1 
O  | 

o 

*% 

«.° 

az  0 

_ 0_ 

Sr  0 

*«r' 

dcl° 

_ 

d  Dy 

dx  0 
_ 0^ 

dz  ° 

_ 

ar  0 

do0 

m 

0 

5pm 

$0  0 
m 

For  example 


xk 

Z.c 

i 

,  0 

2pk 

0 

.1 

=  N 

X1 

Z1 

0 

+D 

,  2  2 
dl  *pl% 

d  2-o  0 

Z1 

0 

h 

0 

X 

Z 

0 

l 

0 

c 

?xo 

m 

rc 

P 

e 

m 

9pk  -  D 

5.6  3  Residuals 

Using  approximations  xo°,  z°  rather  than  true  values  will 
cause  the  condition  equations  to  yield  non-zero  values  Therefore,  we 
shall  define  as  the  residual  of  the  i1"*1  condition  function 

16 

g.(x°,  z°)  =  )  a..  U.°-(3.  =  e.  i=l,  .,.14 
i  o  ’  o  ^  rj  j  Fi  r 


U.  is  evaluated  atx=x°z=z° 
j  oo 


5  6.4  Least  Squares  Adjustment 

The  observed  values  M^tk  =1,  ..  56}and  the  2  approximate 
parameters  x^,  z will  be  adjusted  to  satisfy  the  condition  equations 
while  simultaneously  minimizing  the  sum  of  the  squares  of  the  M  adjustments 
5^;  each  observation  is  given  equal  weight  since  the  errors  in  the 
were  assumed  of  equal  variance  (see  Sec.  7.2  ) 
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Consider  the  14  condition  equations,  r4  ,  i=l , . ,  14  (see 
Section  5 .6 . 1 ) 

Let  the  column  vector  of  the  residuals  he  E 


'€1  " 

a  :  As 

_v 


3Si  98i  ssi  ? 

smT  :  B=  I  aiT  FT*  :  b=  . : 

k  °  °  8W 


=  V  be  the  correction  vector  of 


The  condition  equations  can  be  expanded  as  Taylor  series  and 
retaining  only  the  first  powers  and  constant  terms,  they  become 


i,  =  )  a . ,  &,  +  )  b.  v  +  e.= 

1  L  iq  q  t 


or  in  matrix  form 


A  6  +  BV  +  E  =  0 


Define  the  sum  of  squares  to  be  minimized  as  s  and  subtracting 
14  zero  terms  as  follows 


■  s  i  ***  -2  i  v  -i 


M  k-1 


where  the  are  a  set  of  14  Lagrange  multipliers  whose 
column  matrix  is 

To  obtain  a  minimum,  differentiate  s  with  respect  to  5^  and 

and  equate  the  resulting  expressions  to  zero. 

,  .  '  14 


1  AS  J_ 

2  „  2  °k  ■ 

M 


a.. 
i  ik 


1  \  x  , 

2  <>'■>  ~  )  i  iq 
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the  vaiuc  >7 


Let  7  be  a  diagon:)  matrix  each  element  c£  which  has 
2 


These  expressions  may  be  written  in  matrix  form  as 

(o°i  1  ■  aF  \  =  0 

BT  *  =  0 


o  1  T 

*i  =  1  a  ■  A  \  and 

A  b  -  BV  *  E  *  0  becomes 
c  ■  ’  T 

A  ’  >  *  A  >,  ■*-BV  +  E=  0or 
•1  T  , 

».  =  -An  A  (BV  *  Ei  and 
r  o'*  T  1 

a  “;Aci  A  .  tpv  *  E»  =  0 


therefore 


(B1  'Ao°  AT-  1  B.'1  (Br  (An0  A1)"1 


I)  =  B1  (Acs°  ATi  1 


M  ‘  11  B 


(N  1  may  be  a  desired  output/ 


V  =  N  i!  E  is  determined  as  the  column  matrix  of  the 


adjustments  and  the  adjusted  vaiues  of  x,  z  are  obtained  from 


V  -  V  u.  U 


-•  VctSti'y  De*.ermin.ar  lan 

5.7  1  Definitions  and  Assumptions 

I’eiocitv,  as  us’d  he.  e,  is  defined  as  the  average  velociry 
in  the  interval  /^t  immediately  following  passage  through  the  beam's 
central  plane 

In  the  determination  of  velocity,  it  is  assumed  that  the 
adjusted  values  of  the  position  coordinates  are  available 
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573  Residuals 
Define  t 


j  *  14 


j  =  1,2, 3, 4  as  follows. 


&Pi 


M8j  _ 

At  *  Xo  fjx„ 


8Pi 


o 


where  p  =  j  +  14 


Using  approximate  values  x  0  and  z  °  instead  of  true  values 

o  o 

will  cause  the  €.  ,,  to  be  non-zero. 

J  +  14 


574  Least  Square  Adjustment 

Let  v.  ,  v  denote  the  adjustments  to  x  °,  z  °,  and  x  ,  z 
x  z  J  o  o  o  o 

the  adjusted  values  It  is  assumed  that  the  adjusted  values  x  and  z 

J  o  o 

used  are  error- less 

To  find  the  corrections  to  x  0  and  z  0  we  evaluate 

o  o 
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SECTION  6 

IE3I  CASE  FOR  DATA  REDUCTION  MODEL 


6  1  Feasibility  Study 

In  crder  tc  shew  the  feasibility  cf  the  data  reduction  method  de- 

jcribvJ  in  s.'tien  3  ?  typical  case  was  selected  for  simulation  on  a 

computer  This  work  was  in  progress  at  the  termination  cf  the  contract 

The  following  is  »  description  of  the  proposed  test  case,  (sea  flow  chart, 

Figure  3)  At  the  time  cf  cessation  of  work,  the  matrix  involved  in 

p-sition  finding  bed  net  been  inverted  It  cresencly  developed  that  this 

was  due  tc  calculation  cf  errcnecus  coefficients  A  delay  caused  oy  the 

unavailability  - f  the  computer  and  subseouenc  funds  shortages  did  net 

permit  ch.  simulation  cf  the  cypicsl  -ore  to  be  resumed.  As  mentioned  in 

Section  3,  it  is  iicced  chat  :n  extension  in  scope  cf  the  present  work  may 

permit  us  cc  pursue  this  investigation 

Synthetic  data  was  generated  f-r  the  minimum  beam  width  of  8°  and 

x  1 400 ,  :  -450,  0->5c,  and  a  circular  erbit  using  the  methods  cf  Aocendix 
c  c 

B  As  m  an  actual  eese,  an  approximate  position  was  to  have  been  found, 
then  adjusted  Sec  section  5  for  technique 

In  the  calculation  of  Che  position  variance -covariance  matrix  for 
the  least  squares  adjustment  cf  position  in  the  test  case,  the  approximate 
values  c£  x.  ana  z^  "re  used  rather  than,  the  true  grid  point  coordinates 
This  is  dc-e  to  simulate  most  clcsely  operational  conditions  For  velocity, 
we  proceed  as  fellows  Using  the  corrected  values  cf  xr  and  z ,  *o°  and 
zQ"  are  found,  a^d  these  are  then  corrected  using  least  sauares  as  before 
(See  section  5) 
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SECTION  7 

ERROR  PROPAGATION  STUD'. 

7  1  Gcnerai  Discussion 

This  section  uttscribes  the  method  used  In  the  preparation  of  the 
.’/<j  curve?.  i  There  i*.  therefore,  some  duplication  between  this  section 
and  the  discussion  oi  least  square-  ir  Section  5,  Data  Reduction  j  l’o 
permit  the  selection  ot  an  optimum  beamwidrh,  the  position  error  matrix, 

N  and  the  velocity  error  matrix.  S  were  calculated  for  bearawidths 
of  8°  and  10*  3C° 

7  1  1  Results  ^ 

Charts  ot"  -  vs  beamwidth  as  abcissa  are 

N>  M  M  M 

presented  for  contours  of  constant  altitude,  zq,  for  various  combinations 

of  horizontal  di -r lacemcnt  x  and  circular  orbit  inclination  (Figures 

o 

6  to  li  I 

Having  selected  a  beamwidth,  line  charts  showing  the  behavior 

o  * 

of  —  —  .  are  presented  as  functions  of  altitude  (z  ),  holding  x 

°M  °M 

and  0  constant  (see  Figures  to  20.  Charts  are  also  presented  showing 

c  n  ~ 

the  behavior  of  —  ,  —  ,  --  .  —  .as  a  function  of  horizontal  displace- 

i*  f  r r  n 

to  *M  M  M 

ment  (x^),  for  constant  altitude  :z0'-  and  for  orbit  plane  angle  <0; 

(bee  Appendix  C  tor  complete  list  of  Figures  » 

7  1  2  Derivation 

The  variance-covariance  matrix  among  the  set  of  variables 
M^,  is  a  symmetrical  square  array  whose  elements  are  the  expected 

values  of  the  cross  products  of  the  errors  in  the  Ms 
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°M . M .  =  E 

1  1 


[(AM,)2]- 


cw.  is  the  variance*  The  notation 
Ml 


a  ' T  arises  from  considering  the  column  matrix  of  the  M. ;  thus 

MM  1 


"l]  [«1  . S>]  ■  K“i 


M56M1 . .  M56M56 


Suppose  chat  y^  (  j  =  l,2)are  a  set  of  variables  which 

are  functions  of  M, 
k 

yl  =  xo  (Mj,  ...  M56) 

y2  =  zQ  (Mr  ...  M56) 

then  56 

c-1  ay. 

Ay.  -  )  •  ' . ,  A  M. 

j  L  aMi  1 

i=l  1 


M  M 


l  l  % 


i=i  q=l 

expected  values 

M  M 

r"*  r-i  syf  dy 


ayi  9yk 

m'  XT  A  Mi  A  Mq  and  C3kin8 

f  q 


o  _  V  V  Ill  lili  B 

yjyk  ‘  L  L  aMq  \\  ’ 

i=l  q=l  1  M 

2n  rayi  i  rayi  i T 

aYYT  =  Jj3M.  J  CMMT  =  5M. 


Let  the  x  and  z  coordinates  of  the  grid  point  be  denoted  by 

xq  and  zq,  respectively.  ALL  FUNCTIONS  OF  x  AND  z  IN  THIS  SECTION  WILL 

BE  TREATED  AS  THOUGH  EVALUATED  ATx=x)z=z,x=x,z=z.  Note 

o’  o’  o’  o 

for  actual  solution  (including  the  test  case)  :  to  adjust  the  approximate 
values,  x^0,  z°  and  xq0,  z^0,  by  the  use  of  least  squares  methods,  the 

variance -covariance  matrix  is  calculated  for  these  approximate  values. 


PH3LCO 


GOVERNMENT  a  INDUSTRIAL  GROUP 
WESTERN  DEVELOPMENT  LABORATORIES 


WDL-TR1278 


The  values  ol  the  parameters  for  which  0/0  matrices  are 

calculated  as 


x  =  50(150)650  statute  miles 
o 


z^  =  150(150)1050  stacuto  miles 
0  =  0(15)90° 

beamwidth  =  b  w  =  8°,  10(5)30° 

(For  0  =  0°,the  first  run  showed  values  of  0/0  which  are 
probably  the  square  of  the  true  values  ^ '  This  is  due  to  the  symmetry 
of  this  case  which  causes  readings  on  the  two  sides  of  the  central  plane 
to  be  eq,.al  This  effectively  reduces  the  amount  of  information  by  one- 
half  and  also  diminishes  the  number  of  condition  equations  by  one-half. 

It  was  intended  that  this  case  be  re-run  with  a  more  irregular  scale,  but 
it  was  not  possible  to  do  so  in  the  time  remaining.) 

Assumptions 

1  Covariances  „  of  the  observed  quantities  M,  ,  ,M_, 

M  M  ^  1  56 

q  p 

are  zero  for  p  /  q 

2  Variances  0  ^  of  Che  observed  quantities  are 

P  2 

equal  for  all  p  and  will  be  denoted  by  0^,  , 

3  One  observation  is  made  for  each  of  the  quantities 

Mi,.  M5&. 

A  Fourteen  condition  equations  are  used  for  position. 


5  Four  condition  equations  are  used  for  velocity- 

6  Second  order  terms  in  the  Taylor  expansion  of  the 
condition  equations  (written  as  residualsmay  be  ignored 


7.  rhe  variance-covariance  matrix  remains  approxi¬ 
mately  constant  in  t  ons  containing  both  the  approximate  and  adjusted 
values . 


(1)  This  hypothesis  is  made  on  the  basis  of  a  study  of  the  trend  in  be¬ 
havior  of  0/0  as  0  varies  from  75°  to  15°,  Extrapolation  seems  to  in¬ 
dicate  a  good  fit  for  0=0  for  the  square  root  of  the  values  for  the 
symmetric  case 
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7 -2  Position  Derivations 

The  condition  equations  are 
.(I*  .(o'* 

A  *  A  =  g.  (M.j,  xo>  zo>  =  e. 


i  =  1,  ...  14 


Let  E  = 


u 

c14 


II 

i 

3"-i 

II  ik 

=  i 

i  =  1 ,  . .  14  k  =  1 ....  56 


N°te’  aik  ^  ®  only  f°r  k  =  j  +  4  (i-1) 


1, 


14 


i  *  1.. 


*ik  2  (Pj  +  Mij>  AjJ'  +  a(,° 


14 

(i) 

AJ2 


where 


j  k  -  4  (i-1)  and  A.^;  is  the  jq1"^  cofactor  (signed  minor)of 
the  determinant 


Let 


Let 


D  = 


2  2 
(  y,  4-  M  '  .H  ‘ 

'"1  ll'  “l 


(P2  *  M.2,2  -d22 
(C3r  M.3)2  -d22 

(C4  ^i6'2  V 


iq 


"j=i 


c. . 
*J 


iPi 

5  x 


aci 


sci 

&0j 


^_pi " 
9  z 


Pi  +  Mu 


°2  +  Mi2 


°3  X  ^i3 


f-A  +  M 


i4 


5ei 


d  x 

52  1 

!fi 

ax 

L  Pp« 

j=i  J 

Sv'i 

i  =  1, 

3Dj 

j  =  i, 

14 

4 
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:i j  =  2  (pj  +  V  Ail  +  Aj2} 

O  _  A(0)  .(0) 

2  Pj  AjX  -  Aj2  , 


D  =  Id 


djl  dJ2 


z  -  Z. 


*  +  _ L  1  + - 

r  +  r.  r  +  r. 
o  j  o  j 


None  that  (r^)2  =  (x)2  +  (z)2 


j  =  1,  ...4 


(r.)2  =  (x  -  X.)2  +  (z  -  Z.)2 


Then  55  =  CD 


Let  a  be  s  diagonal  matrix  each  element  of  which  has  the  value 

2 

a,„, .  In  our  case  we  assume  the  o,~,  =  ow  are  all  equal  and  set  =  1. 
MM  .  MM  M 

T  ,  o'1  T  -1 

Then  let  H  =  B  (A  o  A  ) 

-  BT  (A  AV1 

and  let  N  =  H  B 


then  N 


i  c  o 
■1  _  XX  xz 

a  o 
zx  zz 


7.3  Velocity  Derivations 

The  condition  equations  are 


SPi  ,  3Pi 

-  x  — <-  -  z  — r-i  = 

o  5  x  o  9  z 


j  =  1, 


P  =  j  +  14 
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where 


M0 ,  =  M,  with  U  =  28  +  j 
8j  k 


and  where 


sp.  ap, 

— and  — -4  -  are  evaluated  for  the  adjusted  x  and  z  . 
3  x  3  z  o  o 


Let  A  a  A 


a 

!  rj 

3M 

8j 

r  =  1,  ...  4 
j  =  1,  ...  4 
p  =  r  +  14 


a  .  =  t —  for  r  =  j 

r  j  A  t  J 

=  0  for  r  i-  j 


Let  B  =  B  ,  b 

v  velocity  rq 


9_pl 

3X  C>2 


'"'D 

Then  b  ,  =  - 

pi  3x 

b  - 

p2  3  z 


r  =  1,  ...  4 

q  =  1,  2 

p  =  j  +  14 


=  B 

v 


T  T  -1 

<AV  Av  >  B 


.]  1 


g  0. .  0  . . 

|  XX  X  Z 


ZX  ZZ 


o  -  a. 

X  z 

The  curves  for  —  and  —  are  straight  lines  on  log  log  paper, 

°M  °M 

(see  Figures  6abcd  to  14abcd) ,  when  plotted  vs  bearawidth,  and  so 
o.  a  • 

x  z  n 

—  ,  —  =  k  (b.w  )  It  turns  out  that  in  all  cases  n  =  -1  and  the 


value  of  U  is  10  times  the  sigma  ratio  evaluated  for  b.w.  =  10 
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APPENDIX  A 
NOTATION 


x 

V 


t 

At 


j 

8i 

A 

E(  ) 

Ei 

Hij’Mk,Aij  0C  \ 


T 

a 


horizontal  dimension  coordinate 
transverse  horizontal  dimension  coordinate 
vertical  dimension  coordinate 
x  coordinate  of  grid  point 
z  coordinate  of  grid  point 
time  parameter 

duration  of  single  time  increment  At 
vehicle  path  point  index  (reading  index) 
station  index 

condition  equation  functions 
least  squares  matrix 
expected  value  of  (  ) 
residuals 

range  sum  differences  corresponding  to  the  j1 
reading  at  the  i';l1  station,  k=4(i-l)+j 

receivers 

transmitter 

orbital  plane  inclination  angle 


beamwidth 

range-sum,  i  e  distance  cransmitter-vehicle-ground. 

variance 

covariance 

range  sums  to  grid  point 

x  coordinate  of  the  itr‘  path-point 

y  coordinate  of  the  i*"*1  path-point 
cVi 

z  coordinate  of  the  i  path-point 
adjusted  values  of  unknowns 


A- 1 
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o  0,0  o 

xo  >  zo  -  xo  •  2o 

approximate  values  of  x  ,  z  ,  x  ,  z 

0  0  O  0 

v  v  v  vz 

adjustments  for  x  °  z  °,  x  °,  z  0  to  obtain  x  ,z 

0  0  0  0  0  0 

Qij 

observed  total  cycle  counts. 

Units  lencths:  miles  I'st-amre'* 

time:  seconds 

angles;  degrees 

A-2 
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APPENDIX  B 

GENERATION  OF  SYNTHETIC  ORBIT  DATA 
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B  3  Selection  of  Period  of  Cycle  Court 

Time  of  entry  is  artificially  determined  by  incrementing  time 
until  an  orbit  point  is  feuna  which  lies  inside  the  inner  beam.  (  A 
point  at  which  the  vehicle  is  acquired  by  all  receivers.)  This  will  not, 
except  in  most  unusual  be  the  actual  time  of  entry 

B  3  l  Calculations  of  at  for 0  ~  0 _ (15  )  7  5° 

Calculate 


'R 


■f". 


2a)  *  (z 


'"1 


1/2 


1  2i_x_10*3 
^  R  +_h  3/2 
R 


o  1  .  1 

t  =  -  sin 
a» 


( z_ 


c  *  tan  'Jr 


(R  +  h'  cos  0 


t1  *  t1  3  -A  (in  this  case  A  =  .2  sec) 


Determine  the  smallest  s  such  that 

y  ' t5 i  <  (C*  ' c  1  f tan  V 


where 


y  t)  =  -<R  +  hi  cos  0  sin  cut 

z  t,  =  b'  +  (x  -2a,  sin  0  sin  cot  +  (-b'+  7.  )  cos  cot 
o  o 


let 


Let 


At  =  t  // 

t.  =  (i  -8,  At 
l 

=  (i  7  At 
=  0 


1  <  i  <  7 
8  <  i  <14 


B  3  2  Cal£ulation_of__Ai_for_0_^_9O 
For  each  of  Che  7  secs 


i  =  0 
o 


=  4Q0  miles  0 


90 


z  =  1  SOU 50)  1050  miles 
o 


Calculate  R  h> 


*  (V2a‘‘ 


(z_ 


b')2  I  1/2 


B-2 
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1  24  *  iO 
^  '  R*h  ’  T'T" 

X  ' 

At  y  tM'1  ;-V 

o 

B  4  Orbit  Generation  Equations 
B  4  1  Fosi tion 

Given  "a"  and  a«=u:ning  a  constant  earth  radius  R  to  each  of  the 

five  stations  ,  ✓ 

,  (j  ?  )t/2  ( ...2)1/2 


b  *  (r2  -j)l/2  (s?  -:a<2)l/2 

b  =-(r2  (2?^  j  1,2 

c  =(R2  a2  Y/2-(r2  (2a,2  J  1/2 

See  Figure  5 

Compute  &S  values  i  =  0,  14 

x.  =  Is  ■  b’  +  z&.  sin  0  sin  ^t.  +  fx0‘2a  cos  o>  t. 


y^  =  (R  -  h  cos  0  sin  a,  t. 


z  =  b1  i x  * 2 a  1  sin  0  sin  o>  t.  r  ( -b 1  +  z  )  cos  <o  t. 

1C  10  1 


B  4  2  Velocity 

Compute  the  tve  values 


x  -  u.  •  -b1  +  z  ■  sin  0 
o  o 

z  =  u,  -2a,  sin  0 
c  o 

where  x  *  x(x  ,  z  : 

O  C  O 


B  5  Generation  of  Range  Sums 

Let  o.  denote  range-sum  Comruce  the  56  values 

>  '  l  f  14  i  <  j  < 

f  2  2  2  "l  1/2  f. 

tj  [  i  i  *  J  L  1  J 


'xi 

V 

( z. 

-2  .< 

2  ^  ,  2 
+  y: 


2  I  i/2 
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Compute  the  5  value? 


r  =  [x  2  *•  z  :  1  1,2  (x  .  0 
o  [  o  o  I  o’ 

cj '  ['*.  V2  +  '  V*  ] 


(x  ,  0,  z  >  is  a  grid  ooint 
o  o 


1  <  j  <  4 


B  b  Range  Sum  Pi  i'forencc 
Let 

^ii  =  ?IJ  3oj 

B  o  i  Ir troauction  of  Error  .rto  ‘'true'1  Range-Sum  Differences 
1  n  c-'Sei  to  simulate  realistic  conditions,  errors  in 
Ajj  may  be  introduce  This  is  usually  done  either  by  means  of  fixed 
bias  errors  or  random  error  generation 
!n  anv  case  let 

:,ij  =  *  °rrorij 

I"  the  error  less  case,  ot  course. 

*ij  ^  MiJ 
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Figure  Number 

2a 

2b 

3 

4 


Gabed  Co  14abcd 


15ab  co  20ab 


21ab  co  28ab 


29ab  to  34  ab 


35abcd' 
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Flow  Chare:  Synchecic  Daco  Generation 
Flow  Chare.  Data  Reduction  -  Position 
Flow  Chart-  Data  Reduction  -  Velocity 
Flow  Chart-  Test  Case 
Beam  Configuration 
Station  Configuration 

oooo 

v  ^  X  Z 

A  plot  of  —  vs  Bcamwidth  for  fixed  x  and  0, 

11  H  M  m 

for  constant  contours  of  z. 

o  o 
X  z 

A  plot  of  — ,  — ,  vs  z  for  fixed  x  and  Peamwidth,  for 
M  M 

constant  contours  of  0 

o  o 

A  plot  of  — ,  — ,  vs  x  for  fixed  Bcamwidth  and  0,  for 
!1  M 

constant  contours  of  z 


«  X  Z  X  z 

A  plot  of  ,  —  vs  z  for  fixed  x  and  0?  fcr 

M  M  M  M 

constant  contours  of  Bcamwidth. 

oooo, 

A  plot  of  — ,  — ,  ~~  vs  z  for  fixed  x  and  0,  for 

M  M  M  M 

constant  contours  independent  of  Bcamwidth. 
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STATION  CONFIGURATION 
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